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Abstract 

■ The bow shocks of runaway stars with strong stellar winds of over 2000 km s _1 can serve as particle 
acceleration sites. The conversion from stellar wind luminosity into particle acceleration power has an 
efficiency of the same order of magnitude as those in supernova remnants, based on the radio emission 
from the bow shock region of runaway star BD+43°3654 (Benaglia et al. 2010). If this object exhibits 
typical characteristics, then runaway star systems can contribute a non-negligible fraction of Galactic 
cosmic-ray electrons. To constrain the maximum energy of accelerated particles from measurements of 
possible non-thermal emissions in the X-ray band, Suzaku observed BD+43°3654 in April 2011 with an 
exposure of 99 ks. Because the onboard instruments have a stable and low background level, Suzaku 
detected a possible enhancement over the background of 7.6 ±3.4 cnt arcmin -2 at the bow shock region, 

t— > ■ where the error represents the 3 sigma statistics only. However, the excess is not significant within the 

\ systematic errors of non- X-ray and cosmic-ray backgrounds of the X-ray Imaging Spectrometer, which are 

±6.0 and ±34 cnt arcmin -2 , respectively, and the 3-sigma upper limit in the X-ray luminosity from the 
shock region, which is 1.1 x 10 32 erg s _1 per 41.2 arcmin 2 in the 0.5 to 10 keV band. This result leads 
to three conclusions: (1) a shock-heating process is inefficient on this system; (2) the maximum energy 
of electrons does not exceed ~ 10 TcV, corresponding to a Lorentz factor of less than 10 7 ; and (3) the 
magnetic field in the shock acceleration site might not be as turbulent as those in pulsar wind nebulae and 

■ supernova remnants. 

lO \ Key words: acceleration of particles — X-rays: stars — stars: early-type — stars: individual 

^ (BD+43°3654) 

o; 

(N ■ 

, 1. Introduction ing the non-thermal universe. 

••■ . . _ , ,. , , ., In this paper, we focus on early- type massive stars with 

. 1.1. Runaway stars as particle acceleration sites , ... , ,, nn , _i , . ^r. 

. !_L , velocities greater than 30 km s , that is, runaway OB 

An important goal of astrophysics is to understand stars, as another source of non-thermal emission. The 

j_h ' the non-thermal phenomena in the universe. In our most prominent characteristics of runaway OB stars are 

5^ ■ galaxy, pulsars and X-ray binaries are well-known ob- their high spatial velocities and their very fast stellar 

jects that generate high-energy electrons and possibly winds. The velocities of the stellar winds reach a few 

protons, which emit non-thermal photons (e.g., Sturrock thousand kilometers per second, which is comparable with 

1971; Ruderman & Sutherland 1975; Cheng et al. 1986a; the shock speeds of young SNRs (Bamba et al. 2005 and 

Cheng et al. 1986b; Harding & Gaisscr 1990). Recently, references therein). Thus, in the rest frame of an OB 

sensitive searches of non-thermal X-rays have led to mag- star, two shocks (a bow shock and a wind termination 

netic white dwarfs also being considered as possible non- shock) and a contact discontinuity appear in order to bal- 

thermal emitters (Terada et al. 2008b; Terada et al. ance the ram pressure of the interstellar gas with that of 

2010; Terada & Dotani 2011). In addition to rotating the stellar wind (del Valle & Romero 2012). Bow shock 

magnetospheres, like pulsars, shocks in supernova rem- structures can often be seen in infrared images of runaway 

nants (SNRs) and pulsar wind nebulae (PWNe) are the OB stars (van Buren & McCray 1988; Peri et al 2012). 

most common Galactic cosmic-ray acceleration sites (e.g., Benaglia et al. (2010) have reported radio observations in 

Shklovskii 1953; Ginzburg 1953; Bell 1978; Lagage & the NRAO-VLA Sky Survey program, indicating the pres- 

Cesarsky 1983; Kennel & Coroniti 1984; Koyama et al. encc of high-energy electrons in the bow shock region of 

1995; Rees & Gunn 1974). Searching for new types of runaway star BD+43°3654. 

particle-acceleration sites is a possible key to understand- The object BD+43°3654 is thought to be a blue strag- 
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gler system formed by a close encounter between two tight 
massive binaries (Gvaramadze & Bomans 2008) in the core 
of the Cygnus OB2 association (Comeron et al. 2002 and 
references therein) , which is the most massive OB associa- 
tion in the solar neighborhood, at a distance of d = 1.4 kpc 
(Hanson 2003; Comeron & Pasquali 2007). The spectral 
type of the star is 04If with an age of about 1.6 Myr and 
a stellar mass of 70 ± 15 Mq (Comeron & Pasquali 2007); 
the star is one of the most massive runaway stars known 
today. The heliocentric radial velocity is measured to be 
= -66.2 ±9.4 km s" 1 (Kobulnicky et al. 2010). The 
typical mass-loss rate for 041 stars with a stellar wind ve- 
locity of v sw = 2300 km s _1 isM = 10~ 5 A7© yr~ x (Markova 
et al. 2004; Repolust et al. 2004). The bow shock struc- 
ture was found by IRAS observations in a systematic 60 
/im survey (van Burcn & McCray 1988), and was con- 
firmed by the Midcourse Space experiment in the D and 
E bands (Comeron & Pasquali 2007). Synchrotron emis- 
sion was also found from this region in the radio bands at 
1.42 GHz and 4.86 GHz, as already mentioned (Benaglia 
et al. 2010). 

1.2. Similarity with SNR cases 

We can investigate the runaway system as a shock accel- 
eration site by using its similarities with SNRs: velocity 
of stellar wind, system size, environment of inter-stellar 
matter, etc. If we adopt the typical parameters of stellar 
winds from 041 stars (Markova et al. 2004; Repolust et 
al. 2004), the stellar wind has a kinetic luminosity up to 

r erg s . (1) 

2300 kms~V 

The termination shock causes particle acceleration which 
generates high-energy protons and electrons (del Valle & 
Romero 2012). In order to account for the radio lumi- 
nosity of BD+43°3654 reported by Benaglia et al. (2010), 
high energy electrons with the energy index of p ~ 2.0 
should have a total energy of E c ~ 2 x 10 45 ergs, under 
the assumption of equipartition between the energy den- 
sities of electrons u c and magnetic field u-q, as is the case 
with many young SNRs (Bamba et al. 2003; Bamba et al. 
2005). In this case, corresponding magnetic-field strength 
B is about a few tens of micro Gauss. Then, the cooling 
timcscalc of electrons emitting synchrotron radiation with 
characteristic frequency v (Reynolds & Kcohane 1999), 

— "'(^"(uk)" 1 ' 2 - < 2 > 

is much longer than the dynamical timescale, 

%~^ = 3x 10 10 s, (3) 

where R ~ 2 pc is the distance to the shock front of the 
bow shock measured by Benaglia et al. (2010). 

Electrons are accelerated and emitting synchrotron ra- 
diation around the shock region. They escape from the 



shock region flowing out along the stream backward of 
the star in a certain time scale. Assuming that the sys- 
tem is in a steady state, the electron production rate L e 
is roughly identical to the loss rate from the system. In 
that case, L c will be described by E c divided by the es- 
caping time scale. Here, we regard the dynamical effects 
should be dominant, since the cooling time scale of the 
electrons is fairly longer than the dynamical time scale, 
Td y , as described above (equations (2) and (3)). Thus 
a rough estimate of L c is ~ E c jT& y — 6 x 10 34 erg s _1 , 
which corresponds to about 0.3% of -B sw given in equa- 
tion (1). Therefore, the conversion from the energies of 
stellar winds into the acceleration power of particles in 
this runaway star system has an efficiency of the same 
order of magnitude as those in SNRs. 

Since about 10% of O-stars exhibit bow shock structures 
(Peri et al 2012), the ~ 2,000 runaway stars in our Galaxy 
can produce a non-negligible proportion of about 10% of 
the cosmic-ray electrons observed at Earth (i.e., 10 38 erg 
s _1 ), if BD+43°3654, with its high-energy electron pro- 
duction rate of 6 x 10 34 erg s _1 , represents the typical 
case of a runaway-star system as a particle acceleration 
site. Therefore, it is important to study physical parame- 
ters of the particle acceleration process in BD+43°3654. 

In the sense of energetics, BD+43°3654 is an SNR- 
cquivalcnt particle acceleration site. Then, the next stud- 
ies should be how much energy the particles gain in the 
system. The maximum energy of accelerated electrons, 
£ max , is measured by the roll-off frequency, v xo \\ 1 of the 
synchrotron radiation. If E max is ~ 8-20 TeV as suggested 
by Benaglia et al. (2010), then, according to equation (2) 
in Reynolds & Keohane (1999), i/ ro ii should be ~10 16 - 
10 20 Hz, which comes into the X-ray band. If the spec- 
trum were simply extended from the radio band so that 
f ro ii = oo, then X-rays should be clearly detected with 
the sensitivities of current missions. In that case the ab- 
sence of X-rays would provide a unique upper limit E max . 
However, high sensitivity to diffuse X-ray emissions is re- 
quired in this case. In this study, we performed an X-ray 
search for possible synchrotron emission from the shock 
region of BD+43°3654 with the X-ray satellite Suzaku 
(Mitsuda et al. 2007). 

The rest of this paper is organized as follows. We de- 
scribe the X-ray observations of this object in Section 2, 
summarize the results of Suzaku analyses in Section 3, 
and discuss our results in Section 4. 

2. Observation and Data Reduction 

2.1. Suzaku observation of BD+43°3654 

Suzaku is the fifth in a series of Japanese X-ray satel- 
lites having sensitivities in the 0.2-600 keV band (Mitsuda 
et al. 2007). It carries two X-ray instruments on board: 
the X-ray Imaging Spectrometer (XIS) for the 0.2-12 keV 
band (Koyama et al. 2007), and the Hard X-ray Detector 
(HXD) for the 12-600 keV band (Takahashi et al. 2007). 
Although both instruments have very low-background sig- 
nals (Tawa et al. 2008; Fukazawa et al. 2009), which en- 
able us to perform high sensitivity surveys of X-rays, the 
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Suzaku XIS is most suitable for our purpose because it is 
capable of imaging and stable backgrounds. 

We observed BD+43°3654 and the bow shock region 
with Suzaku from 2011 April 4 1:53 UT to April 6 2:13 UT 
(OBSID=506004010). To put the bow shock region at the 
XIS nominal pointing position, the aiming point was set 
to (a,<5)[J2000] = (20 h 33 m 40.00 s , +44°03'00.0"), although 
the runaway star is at (a, <5)[J2000] = (20 h 33 m 36.077 s , 
+43°59'07.40"). The XIS was operated in the normal 
clocking mode, without window/burst options but with 
the Space Charge-Injection (SCI) function (Nakajima et 
al. 2008). The XIS consists of three front-side illuminated 
(FI) CCD chips (XISO, XIS2, and XIS3) and one back-side 
illuminated (BI) chip (XIS1). XIS2 was not operational 
due to damage from a meteor oid strike. HXD was op- 
erated in nominal mode: half of the 64 PIN diodes were 
operated with a voltage of 400 V and the others at 500 V, 
and the PMTs were operated at nominal gain. 

2.2. Data reduction 

We used the observation datasets processed by the stan- 
dard Suzaku pipeline version 2.5.16.29, with the calibra- 
tion version (CALDBVER) of hxd20101202, xis20101108, 
xrt20100730, and xrs20060410. We used ftools from the 
HEADAS 6.11 package with XSPEC version 12.7.0. 

The source was detected by the XIS in the 0.5-10 keV 
band at an average count rate of 0.15 and 0.28 cnt s _1 
per sensor in the FI- and BI- CCD cameras, respectively. 
Cleaned events of the XIS data were obtained using the 
standard criteria of the pipeline process. The total expo- 
sure of the XIS was 99.0 ks. 

We did not use the HXD-GSO data because of sen- 
sitivity limitations. According to the INTEGRAL cat- 
alog, there were no contaminating sources in the field of 
view (FOV) of the HXD-PIN detector. The cleaned events 
from HXD-PIN were obtained by using the standard cri- 
teria of the process. The total exposure of HXD-PIN 
was 82.5 ks. The averaged count rate of the cleaned PIN 
data including backgrounds was 0.40 cnt s _1 in the 13- 
70 keV band. The non-X-ray background (NXB) events 
were estimated and provided by the HXD team. We 
used them with METHOD='LCFITDT(bgd_d)' and the 
METHODV='2.0ver0804'. The average count rate of the 
NXB was 0.38 cnt s _1 in the same energy range, and the 
residual count rate is consistent with the level of the cos- 
mic X-ray background (CXB, 0.02 cnt s _1 ) within the sys- 
tematic errors of NXB (Fukazawa et al. 2009). Therefore, 
HXD-PIN detected no significant signals in the 13-70 keV 
band, and thus we do not use the HXD data in the fol- 
lowing analyses. 

3. Analyses and Results 

3.1. X-ray sources in the field of view 

Wc first checked the XIS images in the 0.5-10 keV band. 
The spatial distributions of the non X-ray backgrounds 
(NXBs) on the XIS chips in the 0.5-10 keV band were es- 
timated by the Suzaku ftool xisnxbgen (Tawa et al. 2008) 
using the background database in the CALDB area for 
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the period from 120 days before to 30 days after the ob- 
servation. After the subtraction of the NXB image from 
the sky image with the XIS in the same energy bandpass, 
the vignetting effects of the X-ray telescopes (XRTs) were 
corrected by a simulated map calculated for a uniform sky 
input with the ray-tracing tool xissim (Ishisaki 2007). In 
this simulation, we assume that the energy spectrum of 
incident photons follows the power-law shape with a pho- 
ton index of 1.0. Changes of ±1.0 in the photon index 
correspond to changes in image shapes at the ±5% level. 
The vignetting-corrected X-ray image from the XIS in the 
0.5-10 keV band is shown in Fig. 1. 
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Fig. 1. The color map shows the X-ray image of 
BD+43°3654 with the Suzaku XIS in the 0.5-10 keV band. 
The areas irradiated by the calibration sources of 55 Fe are 
excluded from the image. The radio image in the 1.42 GHz 
band obtained by the VLA observation (Benaglia et al. 2010) 
is overlaid as the contour image in cyan; the radio data is valid 
in the cyan box. The black circles represent point sources de- 
tected by XMM-Newton, which are listed in Table 1. The po- 
sition and runaway direction of BD+43°3654 are also shown 
in green with the cross and arrow, respectively. 

It is important in diffuse X-ray observations to extract 
(or to estimate the flux of) dim point sources in the field 
of view, to reduce systematic errors caused by fluctua- 
tions of unresolved sources in CXB. The brightest source 
is BD±43°3654 itself; the X-ray image around the ob- 
ject was consistent with the point-spread function of the 
XRTs. In the X-ray light curve, after barycentric correc- 
tion by aebarycen (Terada et al. 2008a), we found nei- 
ther significant flares nor coherent-periodic signals in the 
16 to 10,000 s range. The X-ray spectrum was well de- 
scribed by a MEKAL model (Mewe et al. 1985; Liedahl 
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et al. 1995; Kaastra et al. 1996) with a temperature of 
O.60t5xog keV and a metal abundance of 0.32^o;i7 solar, 
where the photoabsorption of the hydrogen column den- 
sity was iVn = 1.45+o'iy x 10 22 cm 2 . The emission mea- 
sure and X-ray luminosity in the 0.5-10 keV band were 
EM = 6.5±|J x 10 55 cm" 3 and 3.7±^J x 10 31 erg s -1 , re- 
spectively. The value of 7Vh is consistent with the galac- 
tic neutral hydrogen map by Kalberla et al. (2005), and 
other parameters were well understood as being typical 
of X-ray emissions from O-type stars (Naze et al. 2011), 
following the rough relation between kT (keV) and EM 
(cm -3 ) given by log EM ~ 53.9 - logfcT. 

For the further rejection of contaminating sources, we 
checked the archive data of BD+43°3654 from XMM- 
Newton (OBSID = 0653690101, PI= Zabalza Victor, 
starting from 2010 May 8 08:04 UT, 46.7 ks exposure). 
Several faint point sources that contaminate the FOV of 
the XIS arc listed in table 1 and are plotted as black cir- 
cles in Fig. 1. Therefore, in the following analyses, we can 
exclude contamination by point sources with X-ray fluxes 
above S c = 3.7 x 10~ 14 erg s" 1 cm -2 in the 0.5-10 keV 
band. 

Table 1. Point sources detected with EPIC during the 
XMM-Newton observation of OBSID= 0653690101 (ex- 



BD+43°3654). 


The X-ray 


flux is in units 


erg s — 1 cm -2 


in the 0.5 


to 10 keV band. 


RA 


DEC 


flux 


308.5455 


44.047592 


1.2 


308.54706 


44.058725 


0.47 


308.23159 


44.051389 


0.77 


308.5781 


43.947109 


0.63 


308.22298 


44.067196 


0.70 


308.3077 


44.097479 


0.83 


308.27586 


44.050153 


0.53 


308.37145 


43.969254 


0.40 


308.3562 


44.156242 


0.67 


308.46396 


44.016647 


0.37 



3.2. Upper limit of the X-ray flux from the bow shock 
region 

To estimate the possible diffuse-X-ray emissions around 
the bow shock region numerically, we first defined 
three regions: 'BD43', 'BSK', and 'BCD', representing 
BD+43°3654, its bow shock, and the background areas, 
respectively The definitions of these areas are shown in 
Fig. 2. The circular areas around the XMM sources are 
excluded from the BSK and BGD regions; the radii of 
these excluded regions were defined by their luminosities. 
The geometrical areas of BD43, BSK, and BGD were 11.5, 
41.2, and 44.8 arcmin , respectively. Then, we counted 
X-ray events detected with the XIS, taking into account 
vignetting effects of the XRTs, we obtained 288.1 ± 11.9, 
81.6 ±2.6, and 74.0 ±2.3 cnt arcmin^ 2 , for BD43, BSK, 
and BGD, respectively, where the errors are the 99% sta- 
tistical ones. Thus, the events in the BSK region exceed 
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those in the BGD region by 7.6 ±3.4 cnt arcmin -2 statis- 
tically. 
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Fig. 2. The energy resolved X-ray images from the Suzaku 
XIS. The 0.5—2 keV and 2-10 keV images are shown in red 
and blue, respectively. Point sources from XMM-Newton are 
plotted as black circles, and the radio contour in the 1.42 
GHz band (Bcnaglia et al. 2010) is displayed in cyan. The 
definitions of regions, BGD, BSK, and BD43, are outlined in 
green, yellow, and white, respectively. 

The measurement of the enhanced counts in the BSK 
region ( 7.6 ±3.4 cnt arcmin -2 ) has systematic errors in 
estimations of the NXB and CXB. That of NXB is re- 
ported as 5% by Tawa et al. (2008), and we adopt this 
value in this paper. The systematic error in the CXB 
level is caused by fluctuations in numbers of unresolved 
point sources (which are mainly active galactic nuclei). 
According to the systematic studies of CXB with ASCA 
by Kushino et al. (2002), the CXB flux within the field 
of view of the GIS (Gas Imaging Spectrometer; Ohashi 
et al. 1996; 0.5 deg 2 ) fluctuated by 6.5% in the 2 to 10 
keV band. Since the size of the BSK region is about 42 
arcmin 2 , the fluctuation would be 43% for the BSK re- 
gion by the XIS. As described in session 3.1, we removed 
dimmer point sources from the XIS data than the GIS ob- 
servation; i.e., the flux limit for the GIS was S c — 2 x 10 -13 
erg s -1 cm -2 (Kushino ct al. 2002), whereas our analy- 
ses achieved 4 x 10 -14 erg s -1 cm -2 (section 3.1). If we 
adopt the index of the log N - log S relation for the pop- 
ulation of point sources in CXB is 1.5, the CXB flux in 
our observation should have 12% fluctuation after removal 
of point sources with XMM-Newton. Therefore, the sys- 
tematic errors of NXB and CXB are 5% and 12% for this 
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FOV, respectively. Numerically, the systematic errors for 
the NXB and CXB in the count rates of BD43, BSK, and 
BCD were ±1.89 ±20.3, ±1.91 ±19.2, and ±5.69 ±27.8 
cnt arcmin -2 , respectively. Thus, the enhancement was 
7.6 ± 3.4 ± 6.0 ± 34 cnt arcmin -2 , where the first, second, 
third errors represent statistics, systematics in the NXB 
estimation, and systematics in the CXB estimation, re- 
spectively; i.e., the X-ray count rate in the BSK region 
is consistent with that in the BGD region within the sys- 
tematic errors. 

The X-ray spectrum of the BSK region, in which 
the vignctting-corrected CXB and the NXB were sub- 
tracted, can be well reproduced by a power-law model 
with a photon index of l-lO^Q ^glQ Qg and an X-ray flux 
of 2.56t°;? 3 +°™ x 10" 13 erg s" 1 cm" 2 in the 0.5-10 kcV 
band, where the first and second errors, respectively, cor- 
respond to the statistical error and the combined system- 
atic errors for CXB and NXB, at the 3 sigma level. Thus, 
the 3-sigma upper limit of the X-ray luminosity from the 
BSK region in the 0.5-10 kcV band is 1.1 x 10 32 erg s" 1 
at a distance of 1.4 kpc. 

4. Discussion 

We searched for diffuse X-ray emissions from the shock 
region of BD±43°3654 with Suzaku (see Section 2), and 
found an upper limit for the diffuse X-ray luminosity of 
1.1 x 10 32 erg s" 1 per 41.2 arcmin 2 in the 0.5-10 keV 
band (see Section 3.2). By analogy with SNR systems 
(see Section 1.2), if the termination shock occurs in the 
region, interstellar materials would be shock-heated, and 
electrons or protons would be shock-accelerated. We dis- 
cuss the efficiency of the heating by shock in Section 4.1, 
and the acceleration process in Section 4.2. 

4-1- Efficiency of thermalization by the shock 

If we assume that the value of v sw for the object 
is the typical value of 2300 km s _1 (Markova et al. 
2004; Rcpolust et al. 2004) and that the stellar gas has 
a typical temperature of Tj s ~ 10 5 K and acts as an ideal 
gas having a heating ratio of 7 = | , then the sound speed 

in the gas is c s = 4 x 10 6 (y^%) (573) cm s -1 . This 
value is much lower than v sw . Therefore, we can expect a 
strong termination shock to occur inside the contact dis- 
continuity facing the runaway star, while the bow shock 
occurs on the opposite side. In this situation, the cold 
materials of the stellar wind could be shock-heated to a 
temperature of 

wr ah_ 3 ^ m H ( 4v sw A 2 

= H. 2 (-fL-)( ^ ) keV , (4) 

Vo.615/ 1,2300 km s' 1 J V ' 

where jj,, ttih, and k are the mean ratio of numbers of elec- 
trons and nucleons, the mass of hydrogen, and Boltzmann 
constant, respectively. In other words, we can expect 
hot plasmas to exist immediately behind the termination 



shock. 

In the rest frame of the star, the interstellar-matter 
(ISM) gas flows toward the star and forms shocks through 
interaction with the stellar wind. The incoming energy 
of the ISM flow per unit time is (pism^o/2) x (47ri? 2 ), 
where we have neglected enthalpy, and vq is the stel- 
lar bulk speed through the ISM. This gives a smaller 
contribution to the energy source of a possible thermal 
emission than the ejected kinetic luminosity of the stel- 
lar wind, E sw ~ (l/2)Mii 2 w ~ 27r j cw 3 w i?§, by a factor of 
F = v sw /vo ~ 40, assuming vq ~ u*. Here, Pism and p are 
the densities of uniform ISM and the stellar wind just in 
front of the termination shock, respectively, and we sim- 
ply assume the pressure balance, pism^q ~ p^w, an( i a 
spherical shape for the termination shock with radius Rq ■ 
However, the shocked, hot region is compressed to have a 
density F 2 p ~ 1.6 x 10 3 p at most, and thus we can expect 
an anisotropy of the emission from hot plasmas: it should 
be brighter toward the runway direction as in the radio 
synchrotron image (Benaglia et al. 2010). 

The Suzaku observation represents the upper limit of 
the X-ray luminosity at < 1.1 x 10 32 erg s _1 , which corre- 
sponds to the emission measure EM of < 3.8 x 10 54 cm -3 
for plasmas with the temperature of kT = 11.2 keV as ex- 
pected by equation (4). This value of the upper limit of 
EM reduces only by factor 4 even if we change kT by one 
order of magnitude lower. The EM can be described as 
EM = Jn 2 th dV = f? c tijVes 7 ?, where n e ,th is the density of 
thermal electrons, V is the volume of the plasma, n e ,th 
is the averaged value of electron density, Vbs is the vol- 
ume of the bow shock region (1.5 x 10 56 cm 3 ; Benaglia et 
al. 2010), and rj is a density- weighted filling factor of the 
plasma in the bow shock region. If we adopt the density of 
cold electrons n c . c \ ~ 6 cm" 3 (Comeron & Pasquali 2007), 
the upper limit of EM can be written into the following 
constraint; 

fl^y^ojxHr 3 . (5) 

Therefore, the shock-heating process is inefficient in the 
shock region of the runaway star BD±43°3654. If all the 
cold materials are heated into the hot plasma (n c ,th = 
fie,cl)) for example, then a very small fraction of the shock 
region would be shock- heated; i.e., when a top thin re- 
gion emits thermal X-rays as discussed in del Valle & 
Romero (2012), the thickness of the region would be less 
than 2 x 10 15 cm. Otherwise, if the all the shock region 

emit thermal X-rays (n = 1), then only (^^j ~ 3 % of 

cold electrons are heated by shocks. We expect further 
deep survey by a micro-calorimeter onboard the near fu- 
ture X-ray mission, ASTRO-H. 

4-2. Constraints on a possible particle acceleration pro- 
cess by stellar winds 

As indicated by Benaglia et al. (2010), high-energy 
electrons with p ~ 2.0 on average should exist around 
the shocked region and emit synchrotron radio radiation. 
From the X-ray observations obtained with Suzaku, we 
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obtained a 3-sigma upper-limit of the X-ray luminosity at 
1.1 x 10 32 erg s" 1 in the 0.5 to 10 keV band. Fig. 3 shows 
the multi-wavelength spectra from the shocked region of 
BD+43°3654. We have also plotted the energy spectra of 
synchrotron emissions for p = 2.5 or 2.0, assuming that 
the synchrotron emitting region is a sphere with a radius 
of 3 x 10 18 cm and that the energy densities of electron 
and magnetic fields are in equipartition (i.e., u c = Ub)- 
The model calculation code is adopted from Tashiro et al. 
(2009). In the calculation, we assume u c = ub = 6.9 x 10 -11 
erg/cm 3 or 2.7 x 10~ n erg/cm 3 , corresponding to mag- 
netic field strengths of B ~ 42 /tG or 25 /iG for p = 2.5 
or 2.0, respectively. Therefore, to account for the radio 
and X-ray luminosities, the maximum Lorentz factor for 
electrons, 7 maX! should be less than 10 7 , when p is larger 
than 2.0. In other words, the maximum energy of accel- 
erated electrons would be -E max < 10 or 5 TeV for p = 2.5 
or 2.0, respectively, corresponding to a roll-off frequency 
of z/ roU < 10 17 Hz. 




Energy (eV) 

Fig. 3. The wide band energy spectrum of the shocked re- 
gion of BD+43°3654 from the radio to the TeV band. The ob- 
servation data from VLA and Suzaku arc plotted in blue. The 
red and magenta lines show the energy distributions of syn- 
chrotron and inverse- Compton emissions from electrons with 
a Lorentz factor of lto 2.0 X 10 7 or 1 to 1.0 X 10 7 , and energy 
index p = 2.5 or 2.0, respectively. In the calculation, we as- 
sume u c = ub~ 6.9 X 10 -11 erg/cm 3 or 2.7 X 10 -11 erg/cm 3 , 
(B ~ 42 /iG or 25 fiG) for p = 2.5 or 2.0, respectively. CXB 
is only the source of inverse Compton emission in this cal- 
culation. The green lines show the sensitivities of the future 
X-ray mission ASTRO-H (Takahashi et al. 2010 and refer- 
ences therein), the Cherenkov Telescope Array (CTA) North, 
and the gamma-ray satellite Fermi. 

Under the diffusive-shock-acceleration (DSA) theory 
(Bell 1978), particles need to be scattered back to the 
shock front multiple times, so magnetic-field turbulence 
in the upstream region of the shock is crucial for the 
DSA mechanism. Let us consider the time scales for the 
constraint of the magnetic-field turbulence, which is in 
the Bohm limit for PWNe and SNRs (c.f., Shibata et al. 
2003; Bamba et al. 2003; Bamba et al. 2005). The cool- 
ing time scale of electrons with the maximum energy via 
synchrotron emission (~ 2.0 kyr and 4.1 kyr for p = 2.5 
and 2.0, respectively, from equation (2)) is longer than the 
dynamical timescale of ~ 1 kyr (see equation (3)). Thus 



£ max is determined from the balance of the acceleration 
and the dynamical timescales, and we obtain 

max ~ 20 £ c dy 

60 / v s \ 2 
= T V2300 km s" 1 / 

where £ ~ (B/SB) 2 is the ratio of the mean free path and 
the gyroradius of electrons (Skilling 1975; Jokipii 1987; 
Bamba et al. 2003; Nakamura 2010). In the Bohm limit, 
£ becomes 1. Thus, the maximum energy should be ~ 
60/£ or 108/£ TeV, for p=2.5 (B ~ 42 /iG) and p=2.0 
(B ~ 25 /iG), respectively. Observationally, these values 
should be less than 10 and 5 TeV, respectively, and so, 
in either case, we find that £ should be larger than ~ 11. 
Therefore, the magnetic field might not be as turbulent 
as in PWNe and SNRs. 

Non thermal emission from protons, as well as electrons, 
could contribute possible Gamma-ray emission in Fermi 
and CTA band (Benaglia et al. 2010; del Valle k Romero 
2012). The low level of turbulence from our results in the 
X-ray band with the assumption of equipartition between 
u e and u m indicates low acceleration efficiency, not only 
for electrons but also for protons. In this situation, the 
maximum energy of protons must be far below the knee 
energy of 10 15 eV. Furthermore, if the magnetic-field tur- 
bulence is generated by accelerated protons themselves, 
as in the case of SNRs (Lucek & Bell 2000), the fact 
that £ > 11 suggests a lower density of cosmic-ray pro- 
tons in this system, which implies the proton injection 
rate is smaller than that in SNRs. 
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